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Abstract. We investigate the existence and meaning of 
a preferred length in the large-scale distribution of 60-/im 
IRAS sources applying three tests to the 2-dimensional 
distribution on the celestial sphere, involving respectively 
the rms fluctuation and the correlation function of the 
source number density, and the peak number statistics. 
For HDM models we find a best value Aq = 45 — 55 h^^ 
Mpc for the peak of the spectral function P(fc), close to 
the preferred wavelength X — 30 — 40 h^^ Mpc of naive 
single-scale perturbation models. A feature should indeed 
be generated in this range by the anticorrelation region of 
^(r), although it is not detected by recent Fourier trans- 
form analyses of redshift catalogs. Our tests provide ac- 
ceptable fits for CDM models with either HqH = 0.3 — 0.5 
or Qoh = 0.1 — 0.2, in agreement with results in current 
literature. However, the quality of such fits is less satis- 
factory than for HDM models, with the exception of the 
peak number statistics, and we are not able to find a pa- 
rameter interval satisfying all of our tests simultaneously. 
These results may indicate that the cosmic spectrum is 
more complicated than the smooth shapes predicted by 
current theoretical models. 
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1. Introduction 

The anisotropy field in the 2-dimensional distribution 
of IRAS sources on angular scales of 1° — 70° was re- 
cently used to derive a characteristic length scale of about 
50h~^ Mpc (Fabbri & Natale 1993; the Hubble constant 
Hq is here lOO/i km s^^Mpc~^), significantly larger than 
the dipole-dominating scale (Villumsen & Strauss 1987; 
Clowes et al. 1987) and comparable to the distance Dq « 
3Qh~^ Mpc at which the galaxy-galaxy correlation func- 
tion ^(r) becomes negative (Guzzo et al. 1991; see also 
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Davis & Peebles 1983). Torres et al. (1994) have shown 
that the Fourier transform of ^(r), suitably modeled to 
reproduce the data of Guzzo et al., should induce a peak 
in the spectral function P{k) of the inhomogeneities at 
a wavelength Apeak — 1-5-Do; further, they find that best 
fitting on COBE data, involving the angular correlation 
function of the cosmic background radiation, provides 
Apeak = 51 ± 18 Mpc. We are thereby lead to sus- 
pect that the characteristic length in IRAS may be close 
to a peak of P{k), although not necessarily to its abso- 
lute maximum. This interpretation, even in its weakest 
form, seems to be not consistent with several works infer- 
ring P{k) from the 3-dimensional distributions of visible 
and IR galaxies and clusters (Vogeley et al. 1992; Fisher 
et al. 1993; Einasto et al. 1993; Jing & Valdarnini, 1993). 
These works, based on Fourier transforms of redshift-space 
distributions, provided spectra which, although not com- 
pletely consistent with each other, steadily increase up to 
or beyond IQQh^^ Mpc: When some peak is detected, it is 
located at large wavelengths, in the range 180 — 300 
Mpc. However, Bahcall et al. (1993) argue that large-scale 
velocities induced by the inhomogeneity field itself, not 
properly taken into account in the above works, may lead 
to largely overestimate the large-scale portion of the spec- 
trum in real (rather than redshift) space. Corrections are 
strongly model dependent, and for some of the CDM mod- 
els considered by Bahcall et al. they shift the peak wave- 
lengths in the IRAS and CfA spectra down to ~ 100/i~^ 
Mpc. Is is thereby important in our view to exploit alter- 
native techniques that are not affected by the problem of 
large scale velocities. 

In this paper we reconsider the question of the signif- 
icance and reliability of the preferred scale that we iden- 
tified in Version I of the IRAS Point Source Catalog. We 
use a larger sample including 12511 sources derived from 
Version II of the Catalog, and we perform three statistical 
tests based on finite-size bins. These include the shapes of 
the functions (i) AA^rms(7,7), the rms fluctuation in the 
source number density vs. angular scale 7 in the range 
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1° — 72°, (ii) C(7,(t), the autocorrelation function of the 

number density at 2° — 25°, and (iii) 7Vpoak(o', T), the num- 
ber of local maxima in the source distribution vs. Gaussian 
beam dispersion a (in the range 1° — 10°) and peak-height 
threshold T . The experimental data are tested against a 
large number (about 2000) of theoretical models. 

First we model the inhomogeneity field by isotropi- 
cally stochastic perturbations endowed with a single wave- 
length A, and search for the best values of A, the cosmo- 
logical density parameter Oq) and the product hV (the 
biasing factor of IRAS sources times the local velocity). 
Although the present sample allows us detecting some de- 
pendency of A on the chosen angular range (in particular 
using the statistics on A'^poak), nevertheless we can identify 
a preferred range around 30 — 40 Mpc, irrespective of 
Oo which cannot be determined by the present tests. We 
then examine continuous spectra modeling both hot and 
cold matter with the usual requirement that P(fc) oc /c at 
large wavelengths. We find that the HDM model in the 
r^o = 1 cosmology gives excellent fits for all of our tests, 
except for one data point of AiVrms at the largest scale 
71.7°. The HDM peak wavelength Aq is only slightly larger 
than A foimd for single-scale models. Individual tests can 
give best values as high as 60/i~^ Mpc with large error 
bars, but two joint tests combining the peak number data 
with rms fluctuation and correlation function data, respec- 
tively, give Ao = 52.5if h'^ Mpc and 45.0t5,5 ^"^ ^pc. 
Although we cannot derive stringent limits on hV from 
individual HDM tests, when we require to simultaneously 
satisfy all of our tests we get a preferred value of 450 — 550 
km/s (in close agreement with results that we find from 
single-scale models for flo = !)• 

Basically the same picture arises for other sharply- 
peaked spectral forms, for which the only important pa- 
rameter is the location Aq of the peak of P(A:). Thus sim- 
ilar results for Aq are found for simple two-power spectra 
with slopes around —1.5 at large fc. The investigation of 
CDM models on the other hand provides a quite differ- 
ent picture. For individual tests we find best-fit values of 
the shape parameter fio^ in one of the ranges 0.3 — 0.5 
or 0.1 — 0.2. Low-density models arc favoured by several 
studies (e.g., Efstathiou et al. 1992). Now for < 0.5 
the peak wavelength is ^ 100/i~^ Mpc; in fact, single-scale 
models are quite unable to represent so broad spectra as 
those implied by CDM models. However, the quality of 
CDM fits is less satisfactory than for HDM models, with 
the exception of the peak number statistics. In addition, 
no range of Qoh is completely satisfactory for all of our 
tests; this latter fact recalls the results of Scharf and La- 
hav (1993). It is well known that pure CDM models do 
not generally provide excellent fits to available data, a fact 
which motivates studies involving mixed matter (Davis et 
al. 1992, Taylor & Rowan-Roinson 1992) or a cosmological 
constant (Efstathiou et al. 1992). However, the compari- 
son of HDM and CDM fits in our work (Sect. 4) is rather 
impressive. 



The spectral window functions that can be suitably 

defined for A7Vi-ms(7, 7) and C{j,(j) show that such tests 
weight small scales less than other 2-dimensional tests 
from other authors (in particular, Scharf and Lahav 1993). 
This provides support to the idea that the preferred scale 
arising from sharply peaked spectra should correspond to 
something of physical. It is also clear that such a feature 
should be regarded as generated by the behaviour of ^(r) 
around 30h~^ Mpc measured by Guzzo et al. (1991). Since 
we do not directly construct an explicit spectral form from 
data, we do not claim that our tests on the distribution 
of IRAS sources necessarily imply an absolute maximum 
of P{k) at scale Aq, nor can we claim that HDM models 
satisfy the bulk of available astronomical data. However, 
our results make us believe that the shape of P{k) be- 
low 100/i~^ Mpc should not be so simple and featureless 
as generally believed, and that a positive slope up to or 
beyond 180h^^ Mpc is probably incorrect. According to 
recent measurements of the angular correlation function 
in the APM survey the spectrum rises above a power law 
around 30h^^ Mpc and seems to be peaked at 100h~^ Mpc 
(Baugh & Efstathiou 1993), a scale smaller than declared 
by most authors on the basis of rcdshift surveys. In our 
view the problem of the detailed shape of P{k) is not so 
clearly solved and deserves further investigations. 

2. Modelling 

Our models assume that the large-scale inhomogeneities 
are described by a isotropic, Gaussian field of perturba- 
tions. The theory of source-count anisotropics has been de- 
veloped (Fabbri 1988, 1992) for a Friedmann background 
with an arbitrary density parameter CIq, perturbed by a 
cosmic density-wave spectrum P{k) such that 
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= / dkk^P{k). 



(1) 



Here the dimensionless wavenumber k is connected to the 
physical wavelength by 

27rc 



(2) 



with S{no) = |1 - r^oP for Oo 7^ 1 and S{1) = 5. AH the 
quantities that will be evaluated theoretically are derived 
from the squared coefficients of a harmonic expansion of 
the source-count anisotropy field AA^ = {N) 5^£,,„ aimYem- 
More precisely, we consider the expectation values 



(3) 



where we sum up the cosmic-structure contribution 
given by 



a|<, = 47r {2i + 1) j dk 



k'^Pjk) 
{k'^-AKf 



(4) 
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and the noise term aj^g originated by discrete sources 



'■Ne 



(21+ I) 
{N) ■ 



(5) 



In Eq. K denotes the sign of space curvature and 
the functions A^{k) are suitable coefficients introduced 
by Fabbri (1988). We could also define window func- 
tions 'i'e(k) as in Scharf and Lahav (1993) by setting 
Ai{k) oc b{P - 4:K)'<i'i{k), with b the biasing factor of 
IRAS sources. 

The analysis of Af(fc) (Fabbri 1992) shows that the 
source properties affect the calculation through the biasing 
factor, the luminosity function, the number evolution fac- 
tor and, to a very small extent, the spectral index. We shall 
provide explicit results for the exponential-logarithmic lu- 
minosity function of Saunders et al. (1990), which best 
fits data with a number evolution factor (1 -I- z)^-'^ . How- 
ever, using the double-power functions of Lawrence et al. 
(1986) or Villumsen and Strauss (1987) does not change 
the conclusions of Sect. 4 (adopting no source- number 
evolution of course, since they were obtained with this 
assumption). Since in any case the assumed source dis- 
tribution extends farther than probed in redshift surveys 
(Strauss et al. 1990), we checked the influence of a cut- 
off placed at some Zmax- Figure 1 displays some typical 
results for aj, showing that we get no significant differ- 
ences provided z^ax ^ 0.1. A sharp cutoff at redshifts 
^max ^ 0.03 would dccrease the amplitude of high order 
harmonics relative to low order ones. Thus we just assume 
that sources still exist at a distance (along the light cone) 
equal to roughly three times the median depth of PSC. 
We do not regard this as an onerous requirement: In prac- 
tice we only assume that the median depth of the catalog 
is determined by the source luminosity function and the 
lower flux limit. The global properties of the universe are 
involved through fig. Further, it was shown that, because 
of the presence of the discrete source term in Eq. (||) , the 
shape of the harmonic spectrum of anisotropics depends 
both on the shape of P{k) and on the product bV, with 
V the local velocity with respect to the cosmic background 
radiation. Since V is generated by the perturbation field 
itself, it is also a measure of the absolute normalization of 
P{k). (This implies that the shape of aj vs. £ should be 
sufficient to fix its own normalization.) However, V may 
be affected by nonlinear disturbances. On the other hand, 
the absolute normalization of aj is affected by high red- 
shift sources more than the shape of the harmonic spec- 
trum (Fabbri 1988), and it might introduce a systematic 
error on the normalization of P{k). We thereby prefer to 
fit bV and the normalization factor of P{k) as two inde- 
pendent parameters. We then check a posteriori whether 
we obtain consistent results for bag , the biased fluctuation 
in 8h~^ Mpc spheres, computed from these two parame- 
ters when they are both available. For the peak number 
test, only the first method (i.e., derivation from bV) can 
be used. 



Fig. 1. The squared harmonic amplitudes vs. I for single-scale 
models with f2o ~ 1 and hX — 10 and 90 Mpc. Full lines rep- 
resent overlapping curves with Zmax = lO'^, 1 and 0.3. Dashed 
lines refer to Zmax = 0.1, and the dotted one to «max = 0.03 



As to the shape of P{k), here we consider both single- 
scale models and continuous spectra. For the reasons out- 
lined above, single-scale models only depend on three im- 
portant parameters, i.e. A, JIq and bV. For a full investiga- 
tion of their influence we generated a regularly spaced 'lat- 
tice' of 900 models for the adopted optimal set of source 
parameters [i.e., using the luminosity function of Saun- 
ders et al. (1990) and a mean spectral index of —1.395]; 
including also the luminosity functions of Lawrence et 
al. (1986) and Villumsen and Strauss (1987), we investi- 
gated about 1500 single-scale models. Continuous spectra, 
which were investigated only in conjunction with the exp- 
onential-logarithmic luminosity function, were modeled by 
simple analytic forms for hot and cold matter. HDM mod- 
els are described by 



P{k) oc k exp 



(6) 



with kc the cutoff wavenumber. The peak wavelength is 
given by Ao = 2 ■ S^nc [HQS{no)kc] . The number of 
models computed with the spectrum in Eq. ^ was about 
200, and we limited ourselves to the critical density cos- 
mology. We also considered about 200 models with power- 
law behaviours at large k. We analyzed about 200 CDM 
models with 



P(fc) 



1 + 1.7fc + 9.0fc5 + l.or 



(7) 
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Fig. 2. The window functions |*f|^ (full lines) and l*c|^ 
(dashed) for ^lo — 1 and sharp-edge beams, compared to some 
single-harmonic window functions of Scharf and Lahav (1993) 
(dash-dotted) and to a few forms of k^P{k) (dotted lines). 
For \^cf the curves refer to 7 = 2.24° and 17.9°; for 
the curve 7 = 71.7° is added. In both cases the lowest curves 
correspond to largest values of 7 for the (arbitrarily chosen) 
relative normalization. The single-harmonic functions refer to 
/ = 4, 10 and 40, and are computed for the 'optimal weighting 
scheme' with Qq = 0.5, as described by Scharf and Lahav. The 
spectra of HDM models (labelled by II45 and II52) correspond 
to h\o — 45 and 52.5 Mpc respectively, and those of CDM 
models (labelled by C50 and C15) to floh = 0.50 and 0.15 



where k = S'(r2o)fc/(3000 ■ iloh). In this case the rel- 
evant parameter is ^oh, that we let vary in the range 
(0.1 — 1). The corresponding peak wavelength is given by 
hXo ~ 52(Oo/i)~^ Mpc. All of the above continuous spec- 
tra are so adjusted as to reproduce the standard Harrison- 
Zeldovich slope at large scales and exhibit a maximum. 
For sharply peaked spectra, which include HDM (but not 
CDM) models, only the location of this maximum is re- 
ally important. Thus the small-scale behaviour in Eq. (^, 
theoretically predicted by simple modeling of dissipative 
processes in hot matter, can be easily changed without 
affecting our conclusions. (This was tested investigating 
the sharply-peaked two-power spectra.) Although one can 
argue that Eq. (^ cannot reproduce such astronomical 
data as the galaxy correlation function on scales of tens of 
Megaparsec (see e.g. Guzzo et al. 1991), it is the simplest 
one can imagine and turns out to pass our tests quite well. 
Also, for the best fitted parameters arising from the tests 
it gives reasonable values of bag ■ 



For each model we calculated 300 harmonic ampli- 
tudes, in order to provide reliable results down to small 
angular scales. We then derived both the squared fluctu- 
ation in the source number per unit solid angle (for two- 
beam experiments) and the angular correlation function 
smeared over finite beamwidths, which are respectively 
given by the equations 



AiV2„3(7,(T) 



1 



\ " 2 Tp2 



{NY 



2n ^ 



a\Ft{a) [l-P,(cos7)] 



and 



(8) 



(9) 



with 7 the angular separation of beam centers 
(beamthrow), a the effective beamwidth, and Pi the Leg- 
endre polynomials. The coefficients Ff(cr) are affected by 
the beam angular response and are well known for both 
Gaussian and sharp-edge beams. In this paper any theo- 
retical results based on Eqs. (^ and (H) refer to sharp-edge 
circular beams, for which a denotes the angular diameter. 
Equation (^ was used to evaluate AA^rms with cr = 7, 
while using (^ a is kept fixed. In both cases we checked 
that calculations are accurate down to 7 ~ 1°, while at 
7 ~ 0.5° we would need more harmonic amplitudes. 

Although the similarity of Eqs. (H) and (|^) seems to 
imply that these two tests are quite equivalent, this is not 
the case. In the limiting case (7 = 0, the discrete-noise 
terms give a vanishing contribution to C(7, cr) because 

a^|Pf (cos 7) = for 7 7^ 0, and for the reasons ex- 
plained in Fabbri (1992) any dependence on the parame- 
ter hV disappears. For finite values of cr this is no more 
strictly true, but the test based on the shape of C(7, u) 
remains quite insensitive to hV . 

A question which is very important for the main issue 
of our paper is, how length scales are weighted in Eqs. (||) 
and (^). Since ^N'^^J^^, a) and C{j, a) depend hnearly on 
a'j , which in its turn is a linear function of P{k), it is pos- 
sible to define two window functions ^'i? and Vtc such that 
the cosmic-structure contributions to the above quantities 
are expressed as (2/7r) / k^P{k) |*i.,c|^ dfc. [Gf. Eq. (8) of 
Scharf and Lahav (1993).] From Eqs. (|), (|) and (|) we 
get 

i*F(fc,7,f^)r 



27r(fc2 - 4i^)2 
J2i2£+l)A^,F^{a) [1 - P,(cos7)] (10) 



and 



\-^cik,j,a)\ 



47r(fc2 - AKf 

Y,{2e+l)AjFH<j)Pi{cosj). 



(11) 
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Figure 2 reports some results for I'I'fI^ and |5'c|^ in 
critical-density models which illustrate general features. 
(For consistency with beamwidths used in our tests, we set 
(T = 7 and a — 2.24°, respectively, for such calculations.) 
It appears that small scales are weighted in two of our tests 
less than in the work of Scharf and Lahav (1993) using 
the 2-Jy IRAS survey. This is due to the statistical bias of 
the IRAS redshift survey including 2650 objects. Generally 
speaking, we expect a large, flux-limited sample to include 
more far-away sources than a redshift survey. Both |^'_f|^ 
and l^'c'l^ increase monotonically with wavelength. The 
actual contribution of a length scale to the signal depends 
on the spectrum and can be estimated as \'^f,c\'^ P{^)- 
For instance, for rms fluctuations in the angular range 
8° — 35° the signals reach their maxima at A « 35 — 50 
Mpc for a HDM spectrum with Aq = bQh~^ Mpc, and at 
A « 45 - 70 /i"^ Mpc for a CDM spectrum with ^^h = 
0.5. For the same models, the signals fall to 10% of their 
maxima at A « 7 — 17 h^^ Mpc, i.e. around the boundary 
of the linear region, and decrease very rapidly at smaller 
scales. Thus nonlinear disturbances are small. 

We further consider the number of local maxima in the 
density of sources over the sky. We define Npcak{o;T) as 
the number of peaks higher than a threshold T, counted 
in the source distribution after smoothing by a Gaus- 
sian beam with dispersion a. Following Bond and Efstha- 
tiou (1987) we take the threshold height T in units of 
y/C{0, a). The peak statistics is described by 



(27r)5 02 - 



(12) 



e 

^(£-l)^(£+l)(£ + 2)a2^^| 



(13) 



and 



7.^ 



H 2 



H -1 



Y^{e - 1) e ie + i){e + 2) ajF, 



(14) 



where the function / was calculated by Fabbri (1992) from 
the theory of Bond & Efsthatiou (1987), and Fe now per- 
tains to Gaussian beams. Inspection of Eqs. (p^l4) shows 
a highly non-linear dependence on aj , so that it is not 
possible to define a window function. However, the strong 
weight given to higher order harmonics in Eqs. (13) and 
(14) proves that perturbations with small A must have 
a strong bearing on the statistics for small beamwidths, 
so that we expect to be able to provide a lower limit to 



a dominating scale by this test. Taking advantage of 300 
harmonic amplitudes allowed us to evaluate A^peak accu- 
rately down to cr ~ 1°. Since now a denotes a beam dis- 
persion rather than a diameter, the investigation of small 
scales is limited somewhat more strongly than for other 
quantities. 

3. Data Analysis 

We used a galaxy sample based on the IRAS Point Source 
Catalog, Version II, including sources at 60 ^m in the 
flux interval 0.6-20 Jy. The galaxies were selected basically 
using the same criteria as in Strauss et al. (1990), but we 
also set the lower flux limit to 0.6 Jy, the completeness 
limit of the IRAS survey (Yahil et al. 1986), and extended 
the region of avoidance around the galactic plane to ±10°. 
To reduce the residual contamination we also excluded 
other sky regions after source binning; see below. This 
sample differs from the one of Fabbri and Natale (1993), 
not only for the higher number of sources but also for a 
different calibration in Version II of the Catalog, leading 
to smaller spectral indexes of the sources. The average 
spectral index calculated from fluxes at 60 and 25 ^m is 
— 1.44, and the one from fluxes at 60 and 100 fj,m is —1.35. 
For the theoretical modelling in Sect. 2 we adopted the 
mean value —1.395, but the results are quite insensitive 
to this choice. 

Sources were collected in bins spanning equal intervals 
(i.e., in Galactic longitude) and A(sin6) (in the sine 
of Galactic latitude). We considered three resolution levels 
leading to three maps with 128x64, 256x128 and 512x256 
bins, respectively, over the entire celestial sphere. An effec- 
tive angular diameter ab such that (7r/4)cr^ = AZ A(sin6) 
can be defined for the bins. For the above resolution levels 
we have cjb = 2.24°, 1.12° and 0.56°, respectively. For the 
purpose of our analysis clean bins with a vanishing number 
of sources are quite significant, but we must exclude re- 
gions contaminated by galactic emission or not sufficiently 
observed during the survey. Hence to each source map we 
applied a mask excluding bins which are (i) not observed 
at least twice in the survey, or (ii) contained within ±10° 
around the galactic plane, or (iii) confused due to large 
number of sources (such as the Orion complex, Taurus 
and Ophiuchus star forming regions, the Large and Small 
Magellanic Clouds), or (iv) contaminated by cirrus emis- 
sion. Figure 3 gives the map of our full sample of 12511 
sources, obtained by using the above mask. By means of 
this procedure, we constructed 3 matrices with indexes 
{ib, ii) spanning Galactic coordinates, whose elements con- 
tained either the number of sources in the bins per unit 
solid angle N{ib, ii) or a masking flag. These matrices will 
be referred to as 'bin maps' in the following. 

In the flrst part of our analysis we compute all the 
possible differences between the values of N{ib,ii) in ad- 
jacent windows and then the rms value. This amounts to 
calculating the quantity ANi-n-is{l,j) of Eq. (||) for sharp- 
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Fig. 3. Map of the 60 /xm, 0.6 - 20 Jy sample 



Fig. 5. Filled dots give the experimental correlation function 
C(7, 2.24°) for a = at = 2.24°, with error bars computed 
at 1 sigma. The full line describes a hot-matter model with 
hXo = 45 Mpc, Qo = 1, bV = 900 km/s. Other curves refer to 

singlescale models: dashed, hXo = 30 Mpc, flo = 0.9, bV = 700 
km/s; dotted and dash-dotted, same parameters as in Fig. 4 



Fig. 4. Filled dots give the rms fiuctuation AiVrms(7,7) in 
IRAS PSC for at = 2.24° and 1.12°, with error bars com- 
puted at 1 sigma. The curves represent theoretical models with 
parameters selected with best fitting procedures described in 
Sect. 4. The full line refers to a hot-matter model with h\o = 60 
Mpc, no = 1, bV = 500 km/s. The remaining curves de- 
scribe the following single-scale models: dashed, hXo = 40 
Mpc, fio = 0.7, bV = 400 km/s; dotted, hXo = 40 Mpc, 
fio = l,bV = 500 km/s; dash-dotted, hXo = 30 Mpc, Qo = 0.9, 
bV = 400 km/s 



edge beams, except for the fact that observation windows 
are not circular. In order to minimize errors introduced 
by non-circular shapes, we maintain beam self-similarity 
in the generation of larger beamwidths: Starting from the 
above bin sizes we merge adjacent beams so that both 
Al and A(sin6) are doubled at each step. The effective 
beam diameter a covers the range from o-f, to 71.7° at log- 
arithmically spaced intervals. This procedure is similar to 
the one described in Fabbri & Natale (1990). The present 
analysis only differs in the treatment of masked bins in 
the merging process: Composed windows are now rejected 
when more than one bin is masked, otherwise the masked 
bin is filled with the average number of sources evaluated 
on all unmasked bins in the whole map. The results ob- 
tained for (76 = 1.12° and 2.24° are given in Fig. 4 with 
error bars at 1 sigma. Ideally the two data sets should co- 
incide up to a multiplicative factor. The figure shows that 
this is verified to a good accuracy: Quite small deviations 
arise from differences in the masks constructed at different 
resolutions. The data set with (j^ = 0.56° is not reported 
here, since it gives no new results and the first data point 
could not be compared to theory for the reasons explained 
in Sect. 2. 

We also computed the correlation function ata = ab = 
2.24° by means of the equation 
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Fig. 6. The peak number vs. threshold T and Gaussian beam 
dispersion a in IRAS PSC. Data points, normalized to the en- 
tire sky, are represented by mesh nodes. Three distinct meshes 
are given, corresponding to (Tb/ff = 0.93 (a), 0.47 (b), 0.23 (c) 



Fig. 8. The peak number vs. threshold T for uj/cr = 0.47. 
Filled dots represent experimental data extrapolated to the en- 
tire sky, with 1-sigma error bars, at three beamwidths a = 1.2° , 
2.4° and 4.8° (with iVpeak decreasing for increasing cr). Theoret- 
ical curves are given for single-scale models with the following 
parameters: full line, h\ = 20 Mpc, Q,o = 0.9, hV = 500 km/s 
(best fit for this set of data); dashed, h\ = 40 Mpc, 510 = 1, 
bV = 500 km/s; dotted, h\ = 30 Mpc, fio = 0.9, bV = 400 
km/s 



Fig. 7. The peak number for ab/<J = 0.47 from the North and 
South halves of the sample, labelled by N and S respectively. 
Data points, normalized to one half of the sky, are represented 
by mesh nodes 



^ [N{it„ii) - (TV)] [7V(jfc,j;) - (iV)] , (15) 

where the summation is performed over all pairs of beams 
such that the angular separation of center coordinates lies 
in the range (7 — ^(it, 7 + ^crb). The number of such 
pairs is np{^, df,), and 7 = noh with n an integer. Masked 
bins are simply skipped in this computation. Now we can- 
not maintain self-similarity; therefore we excluded two 
Galactic polar cups at |6| > 60° where the stretching of 
bin shapes is strong. Figure 5 reports data in the range 
7 e (2.24°, 24.7°). 

In order to treat the peak number statistics, we should 
take into account that theoretical predictions arc available 
for smooth distributions over the celestial sphere. For any 
distribution of discrete sources we encounter the following 
problem, that moving a sharp-edge beam across the sky 
the detected source number varies in a discontinuous man- 
ner. Therefore we must limit ourselves to smooth beams, 
and in this paper wc shall consider Gaussian beams. In 
principle, for each beam position in the sky we should 
weight all sources in the sample by a smoothing Gaussian 



factor exp 



, with 1? the source angular separa- 



tion from the beam center. In practice it is more efficient 
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Fig. 9. The peak number vs. threshold T for at/u — 0.23. 
Filled dots represent experimental data, with f-sigma error 
bars, at beamwidths a — 2.4°, 4.8° and 9.6° (with A^poak 
decreasing for increasing a). The full line describes the sin- 
gle-scale model with hX = 50 Mpc, ilo = 0.9, bV = 800 km/s 
(best fit for this set of data); the dashed and dotted lines refer 
to the same models as in the preceding figure 

to replace the actual content of each bin by an effective, 
smoothed number of sources according to the formula 





exp 











(16) 



where 1?^^^, = 7 i&oi *;o) is the angular separation 

between cells (ib,ii) and {ibo,im)- Masked bins obviously 
are disregarded. Now a complication arises from the fact 
that the measured number of peaks explicitly depends on 
both a and Cb, while theory assuming a continuous sam- 
pling of the source number-density field only predicts the 
(T-dependence. Therefore, to consistently compare data at 
varying values of a we had to require similar geometries, 
i.e. equal ratios ab/o. Since we considered four logarithmi- 
cally spaced values of a in the range (1.2°, 9.6°), from our 
bin maps we can build three sets of data with Cfc /a — 0.23, 
0.47, and 0.93 respectively. Each set includes a triplet of 
values of ct; for instance, for ajy/a = 0.23 we can collect 
data at 2.4°, 4.8° and 9.6° taking advantage of all of our 
bin maps. 

In practice starting from 3 original bin maps we must 
build 9 smoothed maps, corresponding to specific values 
of a and ab^ and in each of them we search for local max- 
ima of the iVsin distribution. The search was performed 
comparing the source number of each bin with those of 



eight contiguous bins. We handled with the presence of 
holes in the map developing algorithms which extended 
the analysis to larger neighborhoods of potential max- 
ima. The presence of very few uncertain cases does not 
affect our results. The identified peaks were then normal- 
ized to the rms source-number fluctuation calculated for 
each smoothed map, and the distribution iVpoakCo", T) was 
computed for 7 equally spaced values of T between 0.1 
and 1.6. Since the number of peaks should then be extrap- 
olated to the whole sky, we have to evaluate the effective 
area associated to the distribution of detected peaks. We 
adopt two procedures: (a) We only include sky bins con- 
nected to the two big unmasked continents at |6| > 10° 
and surrounded by no more than 5 masked bins, and (b) 
we accept all unmasked bins. 

Figure |6| gives A'^pcaklo', T), extrapolated to the whole 
sky from the entire set of unmasked bins using the first 
method above. Since only data sets with fixed ab/<J are 
homogeneous, the figure shows 3 distinct meshes. For the 
case ab/(T = 0.93 only two values of a are considered be- 
cause data at cr = 0.6° could not be compared with theory. 
Accepting all unmasked bins for the computation of areas 
slightly lowers the peak numbers, but this has no impor- 
tant effect as will be seen in Sect. 4. Figure [t] reports 
-^poak(o', 7") for the case ab/cr = 0.47, derived from the 
North and South halves of the sample. In this figure the 
peak number is extrapolated to one half of the sky. Error 
bars at 1 sigma are evaluated as AA'pcak = v^VpoakT/u) 
with /u the unmasked sky fraction. Such errors are re- 
ported in Figs. H and ||. 

4. Results of Best Fits 

For each of our tests we basically followed the same pro- 
cedure which can be outlined as follows. After selecting 
values for the parameters describing the source properties, 
we are left with two free parameters related to cosmic in- 
homogeneities, i.e. the products h\ (or h\o) and hV . For 
single-scale models one more free parameter is the cosmic 
density parameter Jig; for HDM models we set f^o = 1, and 
for CDM models the spectral properties only depend on 
the shape parameter U,Qh. For each parameter set we first 
minimize with respect to a normalization factor, which 
in the AiVj^j^-^g and C(7, a) tests determines the normaliza- 
tion of P{k). For the peak statistics this factor should not 
be fitted in the ideal case of a pointlike sampling of a con- 
tinuous field, namely for Ub/cF — > 0, but in our case (Tb/cr 
can be as large as 0.93, and we have no warrant that the- 
oretical predictions should reproduce finite-bin data ex- 
actly. However the normalization factor turned out to be 
usually close to unity (within ~10%), and it was always 
so in the neighborhoods of the absolute minima. We then 
obtain 2- or 3-dimensional grids, in which we search for 
X^mn- Statistical errors on each parameter are calculated 
projecting out the extremal points of (Xmin + 1) contours 
on the coordinate axes of two-parameter planes, after one 
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more minimization with respect to the third, miplotted 
parameter when available. The errors quoted in the fol- 
lowing tables are quadratic combinations of such 1-sigma 
errors and grid half spacing. 

Single-scale models 

For single-scale models, adopting the optimal source pa- 
rameters we calculated a 3-dimcnsional lattice of equally 
spaced points in the ranges hX = 10 — 100 Mpc, fio = 
0.1 - 1 and bV = 100 - 900 km/s. Table 1 hsts fitted pa- 
rameter values that we obtained from several tests, as well 
as Xmin compared to the number of degrees of freedom. 
Since errors on Oq often exceeded the range investigated, 
we enclosed the best values within parentheses to declare 
their loss of significance. 

Results reported in the table for the rms fluctuation 
tests (denoted by the letter F) refer to two beamthrow 
ranges, 1.12° - 35.8° and 2.24° - 71.7°, and two values 
of CTb. It was especially important to check a possible de- 
pendence of the best-fitted length scale on the beamthrow 
range. No difference appears in the quoted best values of 
ft,A, but Xmin high in the upper range for the difficulty 
of fitting the data point at 71.7°. All of the rms fluctu- 
ation and correlation function data (see the fourth row 
in the table) are consistent with /lA ~ 30 — 40 Mpc and 
bV ~ 300 — 500 km/s. The parameter bV is correlated to 
fio; if we assume fio = 1 (which is not required by our 
data) we can also derive a preferred range bV ~ 400 — 500 
km/s. Results from the peak statistics reported at rows 
5-10 of the table (denoted by the letter P) are associ- 
ated to a beamwidth value. Each entry is in fact obtained 
using 3 beamwidths (of which the label value is the geo- 
metric mean) referring to data points with a single value 
of (Tb/(7. Thus for instance, the cases labelled by P2.4°(a) 
and P2.4°(b) use each 21 data points with a = 1.2°, 2.4° 
and 4.8° and ab/<T — 0.47. The letters a and b refer to 
the methods for extrapolating the number of peaks to the 
whole sky described in Sect. 3. The test combining all of 
our A'^pcak data is denoted by the full range 1° — 10°. En- 
tries in the 6th and 7th row show that consistent results 
were obtained for the North or South hemispheres (which 
was verified for A'poak tests at all beamwidths) . We do not 
report results for the self-similar case with ab/a = 0.93 
because no significant result could be found after exclu- 
sion of data with a — 0.6°. Inspection of the table shows 
that in all cases we consistently find hX = 20 — 30 Mpc 
and 50 Mpc for smaller and larger beamwidths, respec- 
tively. The former can be regarded as a safe lower limit to 
a preferred scale. 

Figure |l^ gives a few contours in the {hX, bV) plane 
corresponding to Xmin + ^A^: with Ax^ = 2.3 and 6.17, 
which for a Gaussian field should delimit the ellipses at 
confidence levels of 1 and 2 sigmas, respectively [cf. Press 
et al. (1988), Chapt. 14]. The figure shows that contours 
for the P tests do not overlap at 2 sigmas, since the 



Fig. 10. The contours in the (hX, bV) plane corresponding 
to levels 2.30 (full lines) and 6.17 (dashed) above Xmim for ^ few 
tests based on single-scale models. The specific tests involved 
are identified by the curve labels, which are also used in Table 
1 and explained in the text 

preferred length scale show some dependence on angular 
scale. However, all of the other tests (included those not 
reported in the Figure) provide intermediate results. 

Finally, since as already remarked A'peak weights high- 
order harmonics much more strongly than AA'rms(7,7) 
and C{-y,a) do, we treated A^pcak data as quite indepen- 
dent from the other two sets of data. This allows us com- 
bining couples of tests, simply summing up the respec- 
tive x^ fields. Entries labeled as F-fP and C-l-P in the 
last two rows of Table 1 respectively combine data from 
the Pl° - 10° (a) case with data from Fl° - 72° (1°) and 
C2° — 25°. The resulting Xmin ^^"^ only slightly higher than 
the sums of minima in the component fields, which obvi- 
ously refer to different values of the fitted parameters: This 
fact is due to the slow variations of x^ around the min- 
ima. Therefore the rather large values of Xmin i^^ Table 1, 
which emphasize the limits of the single-scale approxima- 
tion, mainly arise from some difficulties that the models 
encounter in each of our basic tests. 

4-2. Continuous spectra 

In the case of continuous-spectrum models after minimiza- 
tion of x^ with respect to the normalization factor we built 
2-dimensional grids in the plane {hXo, bV) or {iloh, bV) 
for HDM and CDM respectively. 

For HDM models the range limits were the same as for 
single-scale models, but the 2-dimensional x^ grid was not 
equally spaced in the hXo direction. (The finest spacing, 
adopted around the x^ minima, was 2.5 Mpc.) Table 2 
reports the results obtained in a subset of our best fitting 
procedures, which is sufficient to illustrate our findings. 
The table shows that the values of x^ are very close to 
the numbers of degrees of freedom for all of the reported 
tests. The best fitted values of bos, derived both from 
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Table 1. Best fit parameters for single-scale models 



Case ' 


h\ (Mpc) 




bV (km/s) 


Xiin/DOF 


Fl° -36°(1°) 




(1.0) 


r.nn+50 
^'-"-'-360 


4.1/2 


Fl° - 72°(1°) 




0.7 ±0.5 


4001^3° 


15.2/3 


F2° - 72° (2°) 




0.7 ±0.5 


400+-° 


14.5/2 


C2° - 25° 




(0.9) 




21.1/7 


P2.4°(a) 




(0.9) 


500±410 


14.2/17 


P2.4°(a, N)§ 




(1.0) 




15.6/17 


P2.4°(a, S)§ 


20tf 


(0.8) 


600 ± 350 


20.8/17 


P2.4°(b) 


20112 


(0.9) 


500tS° 


14.6/17 


P4.8°(a) 




(0.9) 


8001^3° 


24.0/17 


P4.8°(b) 




(0.9) 


800 


24.5/17 


Pl° - 10°(a)t 


30l^° 


(0.8) 




30.0/22 


F+P* 




1.0 ±0.6 


O'-'U_230 


47.3/27 


C+P* 


30tl° 


(0.9) 


4QQ + 310 


51.2/32 



' Label F denotes rms fluctuation tests in the specified beamthrow range; rounded-off bin size is within parentheses. C means 
correlation test in the specified beamthrow range; bin size is 2.24°. P marks peak number tests for triplets of a around quoted 
values. For the meaning of (a) and (b) labels see Sect. 3. 
§ N , S = Northern and Southern half of the sky. 

* Full range of a used. 

* Joint tests combine Pl° - 10°(a) data with Fl° - 72°(1°) and C 2° - 25°, respectively. 



bV and from the normalization factor of P{k), are often 
fairly close to the 'canonical' value 0.7 (Lahav et al. 1990, 
Saunders et al. 1992) although there is some spreading in 
the results. One might have expected to find lower val- 
ues of 6tT8 for HDM models; however, this is not the case 
because our spectra are peaked at low wavelengths. The 
only real limitation for so satisfactory results is the ex- 
clusion of the data point at 71.7° in the case of the rms 
fluctuation test. The behaviour of the HDM model repre- 
sented by the full line in Fig. 4 shows how such a point 
spoils the quality of a fit which would otherwise be very 
good. In general, the quality of fits for any sharply-peaked 
(including single-scale and double-power) models is some- 
what worse at larger angles. Thus in Fig. 5 models do 
not strictly follow the oscillations of C(7,2.24°) around 
7 « 20°, and finally, all of our theoretical distributions 
of density peaks slightly deviate from experiment at our 
largest beamwidth, the measured A'^poak(9.6°, T) being sys- 
tematically lower than predicted for T > 0.8. We could not 
extend the analysis to angles larger than 10° because the 
small number of expected peaks is completely masked by 
Poisson noise. 

For the above reasons, we performed the global test 
F+P for HDM combining the Pl° - 10° (a) and Fl° - 



36° (1°) data. Contours at 1- and 2-sigma confidence lev- 
els for the F+P and C+P tests are reported in Fig. ^ 
These global tests appear to be somewhat complementary, 
setting more stringent limits on bV and /iAq respectively. 
We can satisfy both simultaneously with bV ~ 450 — 550 
km/s and /iAq ~ 45 — 55 Mpc. 

Table 3 collects results for CDM tests corresponding 
to those of Table 2. In all of these tests Vl^h is allowed to 
vary in the range (0.1, 1.2). However in the last column 
we also report Xmin i the minimum chi-square obtained 
in the smaller range (0.5, 1.2). The latter range is suit- 
able to test critical density models, letting h vary; in this 
case the constrained minima are always found at h — 0.5, 
but the fits are unsatisfactory with the exception of peak 
statistics. Allowing for low-density models the situation 
improves substantially: For tests involving peak statistics 
only, Xmin is even slightly smaller than for hot matter. 
However, the chi-square minima identified by our tests 
are spread in the range Vl^h = 0.1 — 0.5 and do not ap- 
pear so mutually consistent as in the case of hot matter. 
The situation is made clear by inspection of Fig. |lj: For 
cold matter the C+P and F+P tests allow us identifying 
an overlap region (yi^h ~ 0.4, bV ~ 600 km/s) only at a 
formal confidence level of 2 sigmas. 
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Table 2. Best fit parameters for IfDM models 



Case 


hXo (Mpc) 


bV (km/s) 


basv ^ 




xLn/DOF 


r 1 — OD (^1 j 

C2° - 25° 
Pl° - 10°(a) 
F+P* 
C+P* 


DU_i4 

52.5tf 


cnn+130 

900 ± 560 
700^}^° 


U.oU 
1.14 
0.89 
0.62 
0.89 


U. ( 

1.23 


z.O/o 

8.3/8 
25.6/23 
30.3/27 
34.1/33 



' Calculated from bV in the preceding column. 

^ Calculated from the normalization factor of P{k). 

* Joint tests combine Pl° - 10°(a) data with Fl° - 36°(1°) and C2° - 25°, respectively. 



Table 3. Best fit parameters for CDM models 



Case noh bV (km/s) Xmin/DOF Xmin/DOF 



Fl° - 


36°(1°) 


o.ioing 




5.7/3 


21.1/3 


C2° - 


25° 


0.50l°;i« 


900l^°o 


13.8/8 


13.8/8 


Pl° - 


10° (a) 


0.35«:}« 




23.5/23 


24.9/23 


F+P* 






600ir 


35.2/27 


45.6/27 


C+P* 






9001"° 


38.2/33 


38.8/33 



^ Starred quantities (last column) refer to best fits limited to the range floh > 0.5. 

* Joint tests combine Pl° - 10°(a) data with Fl° - 36°(1°) and C2° - 25°, respectively. 



In order to decide whether our tests really argue 
against CDM models in low-density cosmologies, we 
should discuss the significance of error bars we attach to 
fitted parameters. The confidence regions derived from 

contours stem from purely statistical errors. (Error 
bars quoted in Tables 1-3 also include grid spacing.) 
As pointed out by White, Krauss and Silk (1993) cos- 
mic variance can be more important than instrumental 
noise in anisotropy experiments such as those on the cos- 
mic background radiation. In our case the 'instrumental' 
noise (i.e., the error in counting extragalactic sources in 
sky bins in a definite flux interval) is negligible in front 
of the noise term AN{ih,ii) — y^N{ii,, ii), which is the 
source of our error bars. Cosmic variance is just this noise 
term. In order to check that this interpretation is correct, 
we attributed the standard cosmic-variance uncertainty 

A{af) = af {I + ^) to all of the models which best fit- 
ted any one of our tests, and derived the corresponding 
uncertainties Acv on the theoretical quantities AN^^^^^, 
C(^,a) and A^poak- We found that our declared errors in 
Figs. 4-^ are typically 0.7 — 1.7 times A^y, and this ra- 



tio is larger as it must when theoretical predictions lie 
below experimental data points. The average ratio '^1.2 

^ fu ^ is consistent with theoretical expectation taking 
into account the unmasked portion of the sky. 

As noted by Scaramella and Vittorio (1993), taking 
cosmic variance into account in the computation of the 
field does not provide exactly the same results as simula- 
tions: Models predicting anisotropics larger than experi- 
ment are affected by a larger cosmic variance, and can be 
rejected less easily than we could infer from inspection of 
the field. From these considerations and the inspection 
of Table 3 and Fig. ^ it follows that low-density CDM 
models are not clearly ruled out by our tests. However, 
very simple HDM models work remarkably better. 

5. Discussion 

The results expounded in Sect. 4 show that a common 
picture arises from both single-scale and HDM models. 
The scale A ~ 30 — 40 Mpc suggested by the former 
is only slightly smaller than the peak wavelength Xq aris- 
ing from the latter (in particular, the value 45 — 55 h^^ 
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Fig. 11. The contours in the {bV, h\o) plane corresponding 
to levels 2.30 and 6.17 above Xmim for two global tests on HDM 
models, respectively labelled by F+P and C+P 



Fig. 12. The x^contours in the (bV, Qoh) plane corresponding 
to the same levels above Xmin in the preceding figure, for 
two global tests on CDM models, respectively labelled by F+P 
and C+P 



Mpc implied by the simultaneous consideration of the 
global tests). The linear component of our local velocity 
(in practice, the velocity of the Local Group) is estimated 
around 500 km/s, in agreement with current views on 
the large-scale velocity field, although not with inferences 
of about 1000 kni/s (cf. Collins et al. 1986). Our basic 
conclusions remain unchanged for different source param- 
eters and sharply-peaked spectral shapes provided only 
that IRAS sources already existed at z ~ 0.1. The quality 
of HDM fits suggests that P(fc) is characterized by some 
length scale smaller than appearing from the analyses of 
redshift surveys. 

Possible objections to this result might involve the re- 
liability of the sample, the ability of our tests to effectively 
probe the full perturbation spectrum and contamination 
from small, non-linear scales. 

As to the first point, we are confident that our sam- 
ple is reliable. Similar results (using single-scale models 



and two of the present tests) were obtained from another 
sample derived from IRAS PSC, Version I with more con- 
servative selection criteria (Fabbri & Natale 1990, 1993). 
Thus there is no chance that our results might be affected 
by star contamination. No problem should arise from the 
finite depth of the sample. Applying the peak statistics 
separately to the South and North halves of the sky we 
reduced the effective linear scale of the sample by a fac- 
tor 23 , but we observed no shift of the best-fitted wave- 
length or the confidence regions to lower scales. Also, our 
models carefully take into account the flux interval of the 
sample, the source luminosity function and the relativis- 
tic expression of the luminosity distance; they explicitly 
allow for the effective depth arising from the interplay of 
such factors. So there is no reason why such a length scale 
should strongly affect our tests, nor why a spectrum like 
(H) should fit data better than equation (|^). [In fact, the 
luminosity function of Saunders et al. (1990) does not even 
allow us identifying a length scale so sharply.] The depth 
of the sample could drastically affect the fittings only if 
it arised from a strong cutoff, totally independent of the 
above factors and effective at redshifts < 0.1. There is no 
evidence for such a cutoff. We previously studied the ef- 
fect of varying the lower flux limit on Version I of PSC 
and found no effects of this kind (Fabbri & Natale 1990). 
We further notice that Jing & Valdarnini (1993) use the 
redshift field of 2-Jy IRAS galaxies to infer a peak wave- 
length of ~ 300h~^ Mpc. The median depth of their sam- 
ple is quite smaller than ours. Thus the different results 
arise from the different techniques, not from peculiarities 
of our sample. 

It is therefore important to notice that the shape of 
the window functions which can be defined for our tests 
do not weight small length scales to a pathological ex- 
tent; actually they weight them less than a 2-dimensional 
analysis based on IRAS redshift survey (Scharf and Lahav 
1993). 

All of the above considerations cannot prove of course 
that our best value for Xq should be the location of the 
absolute maximum of P{k). We are aware that best-fit 
procedures are not equivalent to directly constructing the 
shape of P{k) from data. This consideration is especially 
important because we find that our tests, when applied to 
CDM models, point to larger length scales. This appar- 
ently strange result may also be linked to the fact that 
in 2-dimensional analyses length scales are smeared out; 
under this respect 3-dimensional studies (e.g. Feldman et 
al. 1994) are superior. However, the quality of our HDM 
low-wavelength fits and the arguments discussed in the 
Introduction makes us suspect that current spectra ob- 
tained from redshift catalogs may not represent the cor- 
rect P{k) in real space very accurately. Although there 
is a qualitative agreement between spectra inferred from 
visible and IRAS galaxies and clusters, a close inspection 
of the papers by Fisher et al. (1993), Einasto et al. (1993) 
and Jing & Valdarnini (1993) shows significant differences 
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between the spectra proposed by these authors. A con- 
servative interpretation of our results is that at least a 
spectral feature should exist, not detected in the above 
works. Strictly speaking, a local peak around 45/i~^ Mpc 
appears in Fig. 5 of Jing & Valdarnini, but we believe its 
strength is too weak to justify our findings. A recent work 
by Baugh & Efstathiou (1993) studies the 2-dimensional 
correlation function in the APM survey and recovers P{k) 
by the method of Lucy (1974). Their spectrum rises above 
a power law around 30h~^ Mpc and appears to be peaked 
at 10Qh~^ Mpc. Although the authors find that their re- 
sults may be compatible with previous work, still they 
declare that the discrepancy showed by their Fig. 11 is 
disturbing. 

In our opinion the support for a steady increase of 
P{k) up to very large wavelengths is not overwhelming. 
[Difficulties may arise also in connection with data on the 
anisotropics of the cosmic background radiation (Fabbri 
& Torres 1994)]. The problems arising from the compar- 
ison of results found with different techniques may indi- 
cate that the cosmic spectrum is more complicated than 
the smooth theoretical shapes appearing in the current 
literature. 
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